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Abstract 

— ^ ■ We analyze robust stability, in an input-output sense, of switched stable systems. The primary goal (and 

I contribution) of this paper is to design switching strategies to guarantee that input-output stable systems 

• remain so under switching. We propose two types of supervisors: dwell-time and hysteresis based. While 

OA I our results are stated as tools of analysis they serve a clear purpose in design: to improve performance. In 

pH ■ that respect, we illustrate the utility of our findings by concisely addressing a problem of observer design 

^ I for Lur'e-type systems; in particular, we design a hybrid observer that ensures "fast" convergence with 

"low" overshoots. As a second application of our main results we use hybrid control in the context of 
CNJ I synchronization of chaotic oscillators with the goal of reducing control effort; an originality of the hybrid 

control in this context with respect to other contributions in the area is that it exploits the structure and 
^ chaotic behavior (boundedness of solutions) of Lorenz oscillators. 

q 

^ • I. Introduction 

There exist many good reasons and practical motivations to use a set of controllers for a single plant as 
opposed to one controller: for instance, we may think of a complex system whose dynamic behavior can 
only be described satisfactorily, from a practical viewpoint, by using several models, each corresponding to 
, a mode of the system or, in other words, each model being valid for state values in a specified region of the 
I state space. In such case, it is common practice to implement a group of "local" controllers applied depending 
^ ' on the operating mode. Another common situation is that when design constraints are imposed by certain 
^ . "(sub)optimality" goals: for instance, one may ask for a controlled plant to track an operating point under 
constraints regarding transient performance, speed of convergence (to the desired operating point), robustness 
^ ■ with respect to uncertainties or measurement noise, etc. The term ^wZjoptimality is therefore understood in that 
On . sense: to obtain, for instance, the fastest speed of convergence among a set of possibilities (not necessarily 
all possibilities). 
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In such a scenario, it may result convenient to use a set of controllers, each of which achieves the control 
^ ' goal by respecting one of the constraints. The natural question which arises next is how to trade off the 
■ advantages of each "good" controlled system. For instance, given a controller designed to achieve (the) 
fast(est) convergence and given another one which ensures stabilization with (the) small(est) overshoots, 
how to design a switching rule between the two controllers to achieve both, (relatively) fast convergence at 
a low price (small overshoots). We address this problem via so-called supervisory control -cf. [1], [2] and 
[3, Chapter 6]; that is, the goal is to the design an algorithm to appropriately switch between local given 
controls depending on their properties, thereby combining their characteristics to improve performance and 
possibly to achieve "certain optimality" over larger domains. 

For the case of two controls -one local and another global- the problem of uniting controllers has 
different general solutions and motivations regarding stability and robustness -cf. [4], [5], [6], [7], [8]. 
While theoretically challenging in general, focused problems where switching between a local and a global 
controller brings solutions otherwise difficult or even impossible to achieve are common in the control of 
mechanical systems: see [9], [10] for the problems of set-point and tracking control of robot manipulators 
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under switching, [11] for the problem of "swinging up the acrobot", [12], [13] for a hybrid schemes to 
control chained-form systems; etc. 

On a more general basis, switching among more than two controllers imposes significant challenges to 
analysis since classic stability theory does not apply. For instance, as is well-known, switching between 
stable modes of a closed-loop control system does not necessarily yield a stable behavior -cf. [3]. Over 
the past years there has been an exponentially increasing interest on stability and stabilization of switched 
systems, notably for linear switched systems -cf. [14], [15]. An important trend is that based on Lyapunov- 
like methods as for instance, finding a common Lyapunov function -cf. [16], [17], [18]; invariance principles 
-cf. [19], [20], [21]; geometric methods -cf. [22] to mention a few. 

While it stands as a fact that stable systems do not remain stable under totally arbitrary switching, an 
important aspect of stability theory for switched systems focusses on the analysis and design of the switching 
regime. In the context of supervisory control designing the "switching rule" comes to orchestrate a collection 
of (stable) nonlinear control systems. Among the range of conditions one can think of to impose on the 
switching regime in order to establish stability, one of the most notable is dwell time -cf. [23], [1] and 
its many variants (average, weak, state-dependent, etc.) -cf. [24], [25], [26]. Roughly, the introduction of 
dwell time in the switching rule comes to imposing a minimal amount of time that it is imposed to each 
mode to remain active. Following a similar train of thought, we may mention another method which aims 
at analyzing switched systems as time-varying systems with the property persistency of excitation of the 
switching signal -cf. [27]. 

Implicitly, in the previous paragraphs we have discussed the stability problem for switched systems with 
respect to the trivial state motion x{t) = 0. However, irrespectively of the switching method, the control goal 
may vary: from pure stabilization to robustness with respect to external inputs -cf. [28], [29] or measurement 
noise -cf. [5] to the more general setting of output stabilization -cf. [30], [31]. In such case, the stabilization 
goal is to bring an output motion close to a desired operating point. Moreover, it appears natural to consider 
robustness aspects hence to study stability in an input-output sense. This is the context in which we place 
our main results. 

For a given set of "local" controllers, each achieving a control goal in one dynamic mode, we address 
the problem of designing a supervisor achieving stability in an input-output sense that is close to that of 
[32], while improving performance with respect to each controlled mode. Each closed-loop system (say a 
plant and a local controller) is assumed to have a common output to be controlled and for which certain 
"optimality" is imposed (fast convergence, small overshoots, robustness, etc.). Stability is measured taking 
into account external disturbances, as for instance in [33] yet, with respect to an output motion as opposed 
to the full state. More precisely, our definition of input-output stability is given in terms of ICC estimates and 
a generalized norm for the inputs which includes the usual Euclidean and Lebesgue norms. As an aside, the 
definition of input-to-output state stability that we use covers other definitions previously proposed in the 
literature. Under such setting we show that by appropriate switching (hysteresis-based and dwell-time-based) 
input-to-output stable systems conserve such property. 

The scenario previously described is similar (in spirit) to that studied in the recent paper [34] where the 
primary goal is to improve performance of switched control systems; in this reference it is assumed that 
the switching rule is given and the freedom left to the designer resides on setting adequately the initial 
conditions of the controllers. Another recent article along similar lines is [35] where the authors consider 
the more particular case of switching between two controlled systems, one global backstepping controller 
and one local linear-parameter-varying. 

In contrast to works on optimal control of switched systems -cf. [36] and references therein, we do not 
derive (sub)optimal control laws but rather, we assume that a set of controllers is given and then, we design 
a supervisor to orchestrate the switching among these "local" controls. Freedom is left in the way the output 
motion space is partitioned to assign a controller to each of the identified modes. The latter is done on a 
case-by-case basis, following a rule of thumb which is what is often done in control practice. 
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We illustrate the latter point, and thereby the utility of our theoretical findings, by addressing the problem 
of hybrid observer design for Lur'e-type systems and synchronization of chaotic (Lorenz) oscillators. In 
the first case we design a hybrid dwell-time observer with the "optimization" goal of reducing the time of 
convergence of estimation errors. For synchronization of chaotic oscillators we present a hysteresis-based 
supervisor which aims at reducing control effort, as measured both in absolute amplitude (peaks) and energy 
(integral of squared value). The latter is achieved by applying no control over certain regions of the state 
space i.e., using the fact that Lorenz systems have a complex attractor hence, the solutions tend to a bounded 
(relatively "small") domain without control effort. 

The rest of the paper is organized as follows: Section II contains basic definitions and notation; in Section 
III we lay the standing assumptions for our main results, which are presented in Section IV. In Section V-A 
we present an application of our findings to observer design and synchronization of chaotic oscillators. We 
conclude with some remarks in Section VI. 



II. Preliminaries 

As usual, a continuous function a : M4. — > M+ is of class JC if it is strictly increasing and cr ( 0) = 0; 
additionally it is of class /Coo if it is also radially unbounded; a continuous function f3 : IR+ x M+ — > M+ is 
of class ICC, if f5{-,t) G K, for any t E M+, and /3(s, •) is strictly decreasing to zero for any s G 1R+. 

Consider a family of systems 

x = fi{x,d), y = h{x), (1) 

where x £ denotes the state vector; d G denotes a disturbance; y £ W denotes an output and 
i is an index taking values from the countable set X C Z+. Assume that the functions /j : —>■ M", 

h : R.P are continuous and locally Lipschitz with respect to a;, i G Z. It is assumed that d : IR+ —>■ M™ 

is Lebesgue measurable and essentially bounded, i.e., for all i > 



ess sup { |d( t )|, t G [toi } • 



For t = +00 we write ||d|| = ||d||[o,+ oo)- We denote by Mr™ the set of functions such that ||d|| < + 00. 
Our definition of robust stability is stated in terms of the generalised norm S : Mr™ x ^ M_)_ defined, 
for any to ^ by 

t 

S[d,to,t] := a j uj{\d{T)\)dT + b\\A%a)^ (2) 

to 

o, 6>0, a + 6>0, a;G/C. (3) 

The set of essentially bounded functions d s.t. S [d, 0, + cxo] < + 00 is denoted as Mm™ (Mr™ c Mr™). 

Let i : R+ — > T be piecewise constant and continuous from the right then, the family of systems (1) 
defines the following switched system 

= /i{t)(^>d)' y = h{x). (4) 

We say that the switching signal i{t) has an average dwell-time < r/j < +cxo if between switches, for 
any time instants t2>ti>Q it holds that 

for an integer 1 < iVo < + 00 and where iVjj^^^) is the number of discontinuities (switches) of the signal 
i{t) -cf. [24]. If the interval between any two switches is not less than td then the switching signal is said 
to have dwell-time property and A'^o = 1. The system (4), for signal i{t) with average dwell-time or simple 
dwell-time, has a finite number of switches on any finite-time interval and its solution is continuous and 
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defined at least locally -cf. [3]. The switched system, for a switching signal i{t), is called forward complete' 
if for all initial conditions Xq £ M" and inputs d G Mjr™, the solutions x{t,XQ,d) of the switched system 
(4) are defined for all t > 0. We also denote the outputs by y{t, XQ,d) = h{x{t, xq, d)). On occasions we 
may use the shorthand notation x{t) = x{t,XQ,A), y{t) = y{t, XQ,d). 

Definition 1 We say that, for a fixed i G X a forward complete system (1) is Input-to-Output Stable (lOS) 
with respect to the output y the input d and the norm S if, for all Xq E M" and d E Mr™, there exist 
functions j3i E JCC, 7j E /C such that, for all t > 0, 

\y{t,xo,d)\ <(3i{\xo\,t)+-f,{S[d,0,t]). 

We say that a switched forward complete system (4) with i : R_|„ — > X is lOS with respect to the output y, 
the input d and the norm S if, for all xq E M" and d E M^m there exist functions /? E JCC, 7 E /C such 
that, for all t > 0, 

\y{t,xo,d)\ <f3{\xo\,t)+-f{S[d,0,t]). 

Definition 2 We say that for a fixed i E X the forward complete system (1) is state-independent lOS (SIIOS) 
with respect to the output y, the input d and the norm S if, for all xq E and d E Mr™, there exist 
functions E fCC, 7- E /C such that, for all t>0, 

\y{t,xo,d)\ </3'i{\h{xo)lt)+-,i{S[d,0,t]). 

We say that the switched forward complete system (4) with i : M_|_ — > X is SIIOS with respect to the output 
y, the input d and the norm S if, for all xq E M" and d E M]r™, there exist functions /?' E /C/2,7' E JC 
such that, for all t > 0, 

\y{t,xo,d)\ <p'{\h{xo)\,t) + -f'{S[d,0,t]) 
The systems are exponentially SIIOS if I3[{s,r) = ase~^^ or I3'{s,r) = ase~^^ for some a > 0, 6 > 0. 

For the case when S[d,tQ,t\ = ||d||[4jj j), respectively S is defined by (2), closely connected input- 
output stability properties and relations between them for nonlinear dynamical systems can be found in [32] ; 
Lyapunov characterizations of these properties are presented in [38]. The difference between lOS and SIIOS 
properties defined above strives in the dependence on initial conditions: for an SIIOS system, if the initial 
amplitude y(to) is "small" then the amplitude y{i) during transient remains "small"; for lOS systems, even 
for "small" initial amplitudes |^/(to)| the norm of y{t) may be large during transients for "large" values 
of |a;(to )|- However, the asymptotic behavior of lOS and SIIOS systems is the same as they both "forget" 
about about the initial conditions and converge to the set where the output is zero (in the case that there 
are no external disturbances) or to some neighborhood of the latter (where the size of the neighborhood 
is proportional to the disturbances amplitude). For y = x both properties are reduced to the well known 
properties of input-to-state stability (ISS) and integral ISS properties -cf. [39], [40]. 

III. Problem statement 
Our main results rely on the following hypotheses. 

Assumption 1 For each fixed i E X, the system (1) is forward complete and SIIOS with respect to output 
y and input d with functions /?j E ICC, 7^ := 7 E /C and norm S. 

From Definition 2 it holds that /?j(s,0) > s for all s > 0. Without loss of generality we assume that 

/3i{s,0) > s, s > and denote Xi{s) = P^^{s,0). 

'Necessary and sufficient conditions for a dynamical system (1) to be forward complete with S[d,to,t] = ||d||[t(,,t) can be 
found in [37]. 
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The second standing assumption is that we dispose of an appropriate partition of the real line for normed 
output values and to each partition a stable system is assigned. It should be clear that the number of systems, 
A^, and the number of partitions, M, is in general not the same. Let q G [0, M] define the strictly increasing 
sequence {A^} and let each system from the family (1) be labeled 9q e I for each such q. 

M 

Assumption 2 A partition R+ = |J [Ag, Ag+i), Ag = 0, Am+i = +oo is given and, for each q G 

g=0 

[0,M] there exists 9g £ I such that l3o^{Ag+i,0) < Ag+2- Furthermore, let Tg : [A^, A^+i) M+ be 
given continuous, separated-from-zero (i.e., < Tmin = min | inf iTa(s)}} < + oo), bounded 

functions. 

Assumption 2 holds, for instance, if /3g (s,r) < /3i{s,r) for all r G [0,Tq{s)), s G [A^, A^+i), i G X, which 
implies that the output trajectories of the system 6q, starting off in the interval [Ag,Ag+i), converge the 
fastest to zero as compared to the outputs of any other system of the family (1). In such sense, this introduces 
domains of "local optimality" for each system of (1); therefore, by ensuring the "right" switching between 
the locally optimal systems it is possible to guarantee fast (suboptimal) convergence to zero for the output 
of (1). A graphical illustration of Assumption 2 is presented in Fig. 1. Assume that the set / = { 1, } 

00 = 1 9,=3 02 = 1 eM^i = 2 9m = N 

I \ — \ — I I \-\- — k I \ ^ 



Ao = Ai As A3 A^ A,.+i Ar+2 Am-i Am \y\ 

Fig. 1. Illustrative example when Assumption 2 holds 

with > 1. Then, a possible partition of depending on possible values of \y\ is shown above. For 
each interval [ A^, Ag+i ), < q < M the number 9g € T and the function Tg : [ Ag, A^+i ) — > IR+ are 
given. According to Fig. 1 on the intervals [0, Ai ) and [A2, A3) the system x = fi{x,d) from family 
(1) provides "optimal" performance in terms of output speed of convergence, on the interval [Ai, A2 ) the 
system i = /3( a;, d ) possesses the same property and so on. Dwell-time functions Tg : [ Ag, Ag+i ) —>■ M+, 
< q < M are used to ensure the existence of delays between switches and right-continuity of the switching 
signal i{t). In this way each system G X is active for, at least, Tg units of time. 

Remark 1 Note that one can define functions Tg : [Ag, A^+i) —>■ as solutions of 

/?e„(Ai,ro(s)) = 0.5/9e„(Ai,0), sG [Ao,Ai); 

peAs,Tg{s)) = xe,^A\) 
s G [Ag,Ag+i), q G [0, M]. 

Such function Tg estimates the time that is required for the output trajectories of the system 9g, to converge 
to the interval [Ag„i,Aq) from the interval [Ag,Ag+i); the first function Tq is thus artificially imposed. 
One may replace these functions with any other Tg : [Ag,Ag+i) — > R+, < q < M admitting the same 
requirements of continuity, boundedness and that they are separated from zero. For example, a possible 
choice of Tg for all < g < M and for all s G [Ag, Ag+i) is one which solves: 

(3eAs,Tg{s))+j{Dg) = Ag+,, 

max{ 0, Ag+i - (A^, 0) } < -fiDg) < Ag+i . 

In this case, the functions Tg majorate the amount of time that the output trajectories of system dg remain 
in the interval [Ag, Ag+i) while under the influence of disturbances with amplitude smaller (with respect 
to the norm operator S) than Dg. That allows one to exhibit additional switches caused by disturbances 
smaller in norm than Dg. 



5 



Our control objective is to design a supervisor that generates a piecewise constant switching signal i : 
M_i_ — > T (continuous from the right) such that the closed-loop system is lOS with respect to the output y, 
the input d and norm S. Additionally, it is required to use the system 9q in the domain [Ag, A^+i) for each 
< q < M which is assumed to provide local optimality. We address such problem by designing two types 
of supervisors: dwell-time and hysteresis based. 

IV. Supervisory output control 

A. Dwell-time supervisor 

Informally speaking, the dwell-time supervisor is defined as follows: consider a set of lOS systems labeled 
9q with q E [0, . . . , M] and a given partition of the non-negative real line as described in Assumption 2 and 
illustrated in Fig. 1. After a given instant tj at which the output trajectory lays in the interval [Aj, Aj+i) 
the next switching is determined as the earliest moment, passed a dwell time, such that the output in norm 
belongs to a partition [Afc,Afc+i). Note that the kth and the jth intervals are in no particular order with 
respect to each other and are not necessarily contiguous; e.g., the output trajectory \y{t)\ may pass from 
the interval [Ai,A2) to the interval [AsjAg) or viceversa by 'crossing' intermediate intervals during the 
dwell time. Once the switching instant is identified, the switching rule (the supervisor) takes the value 
corresponding to the interval in which the switch occurs; i.e., the kth mode is selected. To avoid undesirable 
infinitely fast switching, besides the (variable) dwell-time, rather subintervals strictly contained in [Aj, Aj+i) 
are considered. 

More precisely, the dwell-time supervisor is defined as follows: 

to = 0, i{to) = r, r G {0, . . . , M}, \h{x{to))\ G [A,., A^+i) ; (5a) 
4:= arginf {\h{x{t))\ € [Ak,XeA^k+i))} , k e {0, .., M }\{i{tj)} 

tj+i = min {4} ; (5b) 

fcG{0,..,A/}\{ife)} 

i{tj+i) = k such that \h{x{tj+i))\ G [Ak,xed^k+i)) C [A^, Afc+i) . (5c) 

The dwell-time depends on the output trajectories i.e., Ti(^tj) = ^i{tj)(l'*(^(tj))l)' where j = 1,2,3, .. . 
are switching times, j is the index of the last switch and the signal i{t) is piecewise constant for all t 

M-l 

such that \h{x{tj^i))\ G = |J [xe, (A^+i), Ag+i). The definition of the set M prevents from a fast 

g=0 

switching phenomenon (chattering regime) and it plays the role of hysteresis -cf. [41], [26], [24]. The 
second mechanism that prevents the system from chattering is output^ dependent dwell-time T^f^^^y Since 
by construction the functions Tq, < g < M are bounded and separated from zero the system (4), with 
supervisor (5), has dwell-time Tq{s) > Tmin and the system undergoes a finite number of switches on any 
finite-time interval. 

Remark 2 Opposite to the supervisor considered in the state regulation problem -cf. [8], [12], in the present 
context, the hysteresis set M is not sufficient, in general, to guarantee a finite number of switches. This 
is because the output regulation set M may be non compact which stymies the estimation of the system's 
behavior on the set e.g. , it is possible that a system's state trajectory (in norm) escapes to infinity while the 
output remains in M. In section IV-B a hysteresis supervisor is studied in further detail. 

Remark 3 The function Tq for the supervisor (5), can be chosen identically equal to zero. In such situation, 
the system has average dwell-time Tmin with Nq = 2. Indeed, according to (5), i{tj) / i{tj^i) so the system 
^0 may become active each second time only. Therefore, the time interval between two switches is bigger 
than Tmin; in this case, the time Tmin is calculated for < g < M. 

^In this respect we mention [25] where dwell time depends on the whole state as opposed to an output of the system. 
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Remark 4 Using Assumption 1 one can prove lOS stability for the switched system (4) by the applying 
dwell-time or average-dwell time approaches; for this, it is required that the dwell-time Tmin be sufficiently 
large -cf. [42] for details). However, as it was shown in [42] such approach is applicable only for locally 
exponentially stable systems (1). Additionally, taking large values for the average dwell time or dwell-time 
constants implies that all functions j3i from Assumption 1 decrease to zero rather slowly (by construction of 
Tmm) which is not desirable from a practical point of view. Such restriction is opposite to our requirements 
since it is desirable to apply locally optimal controls in their domains. 

Theorem 1 Let Assumptions 1 and 2 hold. Then, the system (4) with supervisor (5), measurable disturbances 
d G M]8m and initial conditions to = 0, a^o G M", 

(i) is forward complete and, for all t > Iq, 

\y{t,xo,d)\ < /9(max{2AM,|^(aJo)|},0) +7(5[d,0,t]) (6a) 
P{s,0) := supft(s,0) (6b) 

7(s) := 7(s) + /3e..(27(s),0); (6c) 

(ii) furthermore, if the system dynamics undergoes a finite number of switches over infinite time i.e., if there 
exists N < oo such that {t^} t]\j < ^ (define tAr+i := +oo) then, there exists a continuous function 
(3^ : M+ X M+ M+ such that (3^ {-,8) is strictly increasing for each s and P^{r, •) is strictly decreasing to 
zero for each r, such that^ 

\y{t,xo,d)\ < p''{\h{xo)\,t) + j{S[d,0,oo]); (7) 

(iii) finally, if d(t) = (with no assumption on the number of switchings) the output trajectories y{t) satisfy 

\yit,xo,d)\<p''{\h{xo)\,t) (8) 

where /3° G ICC. 

Proof of (i). On any time interval [Tg, Tg) with Tg > > a finite number of switches N^t,,t^) with the 
upper estimate 

No + iT,-T,) 

J-^[T,,T,) S Tf, 

^ min 

may occur. Between switches, the system's dynamics is continuous and is equivalent to that of a system 
from the family (1) and which is forward complete for i G / fixed arbitrarily. Since the signal i(t) remains 
constant over [Tg, T^.) the solutions of system (4), (5) are continuous and are defined at least locally on 
the same interval. From continuity of solutions of (4), (5) and forward completeness of i; = /j(a;,d) with 
fixed i, it follows that solutions of (4), (5) are also defined at Tg. The same arguments hold for any interval 
[Ts,Te) therefore the switched system (4), (5) is forward complete. 

We prove next that the bound (6) holds. From (5) we have i(0) = r where r G [0, M] and |fi(a;o)| G 
[Ar, Ay+i). According to Assumption 1, in this case, 

\y{t,XQ,d)\ <Mh{xQ%t)+^{S[d,{)M) VtG[0,ti) (9) 

where ti is the first switching time and, according with (5), ti is when the output trajectory belongs to the 
interval [l\i(^i^-^,xe^ti){^i{ti)+i)) with i{ti) G [0, M] and i{ti) / i(0). If ti = +oo the proof ends replacing 
the bound in (6) with 

\y{t,XQ,d)\ </3e^(|^(a;o)|,0)+7(5[d,0,oo]) Vt>0. 

Otherwise, if ti < +oo, by continuity of solutions, continuity of h and forward completeness the solutions 
may be continued up to ti and 

|y(ti)| < xe,,, JA^(t,)+l) ■ (10) 

_ J _ J 

That is, p is "of class /C/I" with the exception that f3 {0, ■) may be different from zero. 



7 



Reconsidering the solutions with new "initial" condition ti, x{ti) we obtain, from the definition of the 
supervisor (5), 

\y{t,x{ti),d)\ </3e,,,^,(|y(ii)|,i-ti)+7(5[d,ti,i2]) i G [^1,^2). 
Using (10) and the definition of Xe,(ti) expression above, we obtain 

\y{t,x{ti),d)\ < A,i(t^)+i + 7(5[d,ti,t2]), te [h,t2). 

By definition, at time instant t2 the system's output reaches for the interval [Aj(t2)5 X0,(f2) (^j(i2)+i)) ^^^^ 
i{t2) G [0, M] and i{t2) / i{ti) and |y(t2)| < Xe,(t2) (^«(t2)+i)- Following similar arguments as above we 
see that 

\y{t, x{t2),d)\ < A^(t,)+l + 7(^ [d, t2, ts]), t e [t2, ts) 
Repeating the previous steps we obtain, for arbitrary time instants tj for any j > 1 

\y{t,x{tj),d)\ < Ai(j^.)+i +7(5[d,tj,tj+i]), t G [tj,tj+i) (11) 

<X9,<.^,(A,(t^)+i)- (12) 

Now, if at some step k the system's output (in norm) belongs to the last interval [Am, Am+i) with A^f+i = 
+ 00 then, the maximum output amplitude may be estimated by 

\y{t)\ < Pe,A\yiik)\,0)+7{S[d,tk,tk+i]), te[tk,tk+i). (13) 

According with (11) the possible values of y(tfe) are bounded by^ 

\y(.tk)\ < AM + 7(S[d,0,t]), 

Therefore, for all t G [tfc,tfc+i), 

\y{t)\ < Pe,, ( Am + l{S [d, 0, t]), ) + 7(5 [d, 0, t]). 

By induction, the latter holds for any integer k > 1 hence, for all t > ti. Using the "triangle" inequality 
/30„(a + 6,0) < /3e„(2a,0) + /?e„(26,0) and the bound (9) we obtain the estimate (6). 

Proof of (ii). We derive next a bound for the output trajectories which leads to (7). Following the definition 
of the supervisor (5) we have (9). According with (5), ti is the first switching time i.e., when the norm of 
the output trajectory belongs to the interval [Ai{t^^,XB^t^){^i{ti)+i)) with i{ti) G [0,M] and i{ti) / i(0). 
If ti = +00 the proof ends noting that (9) implies (7). 

Otherwise, if ti < +00 by continuity of solutions, of h, forward completeness, and the definition of the 
dwell-time function Tj(ti) -cf. Assumption 2, we have 

\y{ti)\ < pe,X\hM,T^in) +l{S{d,0,h]) . (14) 

Under similar arguments let t2 be the second switching instant, that is the first time such that \y{t)\ G 
[A,(t^),X9,(,„(Ai(j^)+i)) with < i{t2) < M and i{t2) / i{h). We have 

\y{t)\<f3e,,,J\y{h)\,t-ti)+j{S[d,h,t2]) VtG[ti,t2) (15) 

and, if t2 < +00, 

\y{t2)\ < Pe.^,J\y{ti)\,T^ir.) +i{S[d,h,t2]) . (16) 

Using (15) in (14) we obtain 

\yit)\ < /39.„„(/?e.(|^(a;o)|,T„,in) + 7(S[d,0,ti]),t-ti) +7(5[d,ti,t2]) yte[h,t2) (17) 
\y{t2)\ < /3e.„„(/3e„(|^(a;o)|,T^in) + 7(S[d,0,ti]),T„,in) +7(5[d,ti,t2]). (18) 

''where, in view of the definition of the supervisor ~cf. (5b) we exclude the value i{tk) + 1 — M + 1 — 00. 
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Similar arguments hold for any k> I and all t G [tk--i,tk) i.e., 

\yit)\ < Pe,,,^^^,{Pe.,,^_^,o{... o{^^^(\h{xo)\),T,,,i,,)+---+^^^ 

+-f{S[d,tk-i,tk]) (19) 
^t{s) ■■= PeAs,T^in) + l{S[d,0,ti]). (20) 

By assumption there exists N < oo such that {tfc} — > t^r < oo and ^Ar+i := +00; from this and (19) it 
follows that, for all t G [tf^, 00), 

\y{t)\ < P^{\h{xo)\,t-tk)+7{S{d,tk,^]) (21) 

N times 

P''{s,r) := /3,.,,^,(/3,.<,^_^,^{^«(s),r^in)+--- + 7(5[d,0,cx)])),r). (22) 

The function /?''(•, r) is strictly increasing for each fixed r and P^{s, •) is strictly decreasing to zero for each 
fixed s. 

Define k^{s) as the right-hand side of (6a) i.e., n^{s) := /3(max{2AM, «}, 0) + 7(5[d, 0, cx)]). Then, f3'^ 
may be defined as any continuous function satisfying the following: f3^{-,r) is strictly increasing for each 
r > 0; /3''(s, •) is strictly decreasing for each s > 0; 

^^{s,t) > f3^{s,t-tN) yt>tN; 
p'^{s,t) := i3*{s,t)e}CC VtG[0,t7v] 

such that P*{s,tN) = k^{s). 

Proof of (Hi). Finally, consider d = and let the initial conditions be i(0) = r with r G [0, M] where 
|h,(x(0))| G [A,., A^+i). If |^(x(0))| G [xe,X^r+i), ^r+i)^ the output trajectory \y{t)\ may, in general, 
reach the interval [A^+i, xe^+i (Ar+2))- On the other hand, if |/i(x(0))| G [A^, xe^(^r+i)) for any r we 
have, by assumption, \y{t)\ < /36i^(|y(0)|, 0) which implies that \y{t)\ < A^+i i.e., the output trajectory (in 
norm) cannot reach higher intervals than [A,.+i, xe^+i (Ar.+2)) and may only decrease for at least Tq^ units 
of time. Therefore, there exists t' such that the system's output trajectories reach the lower next interval 
[Ar-i, X6»^_i (Ar)). Repeating the reasoning we conclude that \y{t)\ may only continue decreasing to zero 
(right to left on the real line -cf. Figure 1) which implies that the system (1) can undergo only a finite 
number of switches over [0, cxo). The proof ends. ■ 

Thus system (4), (5) has bounded output trajectories for any bounded disturbances and that the system's 
output converges to zero in the case of that there are no disturbances. Opposite to the case where stability 
of the switched system follows for large enough values of dwell-time, here dwell-time Tynin > can be 
chosen arbitrarily small provided that there is a finite number of switches over any finite time interval. 

To take into account lOS systems (cf Definition 1) one needs additional requirements on its detectability 
properties. For example, in [43] it was proved that lOS systems with input-output-to-state stability property 
(a variant of robust detectability property) is input-to-state stable. In terms of the S norm previously defined 
we have the following propositions which follow as corollaries of Theorem 1 . 

Corollary 1 Let all conditions of Theorem 1 be satisfied and let Tg satisfy 

/Je„(Ai,To(s)) = 0.5Ai, sg[Ao,Ai); (23) 

/3,,(s,T,(s)) = xe,_,(A,), sG[Ag,Ag+i), 0<q<M. (24) 

Then, the time Tq.sAi that is required for \y{t)\ to converging into the set {\y\ < 0.5 Ai} for the case 
d{t) = can be estimated as follows (from (5) i{0) = r): 

r 

To.sA, < 5^Tfc(Afc+i). (25) 

k=0 
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If l3eSs,T,{s)) + j{Dg) = xe,_,(Ag), < -f{Dg) < xe,_d\l s G [Ag,A,+i) for < g < M and 
/30„(Ai,To(s)) + 7(Do) = 0.5Ai, < j{Do) < 0.5Ai, s G [Ao,Ai), then estimate (25) holds under 
disturbances d G Mjr™ such that 

G [A„Ag+i), gG [0,M], tG S [d,tl,tl] < (26) 

Proof. It follows using Equations (5) the definition of Tq and invoking Theorem 1. ■ 

The relation (24) defined Tg{-) implicitly; for instance, for Po^ := kse~^ we obtain Tg{s) := ln(A;s) — 
ln(x6» _i (Ag)) while for Tg(s)) oc l/rg(s) we obtain Tg{s) oc s. Estimate (25) provides upper estimation 
on finite time of practical stabilization of the system with respect to set where |y(t)| < 0.5 Ai. 

Corollary 2 Let all conditions of Theorem 1 be satisfied and expressions (23), (24) be valid. Assume that, 
for all q G (0, M] and s G [Aq,Aq+i), there exists t' G [Tg{s),Tq{s) + Tq_i(xe,_i (Ag))) such that 
> (3e,_,{xe,_,{Ag),t-Tg{s)) for all t G [t',Tg{s) +Tg_^{x0,_A^,)}). 

Then, for the case d{t) = the output of the system (4), (5) has the shortest time of convergence to zero 
compared to that of any other system 6g from the family (1), with initial output values in [Ag, A^+i). 

Proof. Conditions of the corollary imply that starting with system q it is necessary to switch on to system 
q — I when the output reaches for X9,_i(Ag) since system ^ — 1 reaches for xSg-alAg-i) faster than the 
gth system. ■ 

Corollary 2 establishes conditions under which the output of system (4), (5) has shorter time of convergence 
than any fixed system from the initial family (1) over its domain of "optimality" [Ag, Ag+i). The overall 
performance, that is for all t is sub-optimal in the sense that it depends on the way local intervals -cf. 
Assumption 2- and functions Tg are chosen; for instance, different stability and convergence properties may 
be obtained. The purpose of the following example is twofold: firstly, to illustrate the utility of Corollary 
2 by choosing the partitions determined by {Ag} and secondly, to provide an insight of the difficulty to 
obtain "true" optimality in the choice of the latter intervals and associated dwell-times. 

Example 1 Consider a switched system of the form (4), (5) with I = {0,1}, i.e., we have only two 
dynamic modes, each of which is stable in the sense of Definition 2 with functions (3q and (3i. The goal is to 
find a set of conditions that lead to an improved rate of convergence -optimality here is understood in that 
sense. We have M = 2 and the partition in Assumption 2 is = [0, Ai) U [ Ai, + oo) -cf. Figure 1. The 
fact leading to optimality (or not) is the appropriate choice of the threshold Ai. Since we are looking for the 
switching rule that minimizes the time of convergence let us introduce the functions 7^ : IR_|_ such, 

that /3fc(s,Tfc(s, A)) = A, s G M+, A G M+, k = 0,1 (it is assumed that Tk{s,A) = for /3fc(s,0) < A) 
i.e., Tfc is the time that takes for the estimate on the output trajectories of system k to reach a given level 
A. For the level of interest here i.e., A = Ai assume first that Ti(s, xo(Ai)) < To{s, xo(Ai)) for s > Ai; 
in such case, it is reasonable to "turn on" dynamic mode 1 as long as the output trajectory remains within 
the upper interval [ Ai, + oo) hence, we set 6i = 1 and the system dynamics is defined by i = fi{x, d) on 
that interval. If on the other hand Tq{s,0.5s) < 7i(s,0.5s) for s < Ai it results reasonable to set = 
hence, the system's dynamics is given hy x = /g(a;,d) as long as \y{t)\ G [0, Ai). After the conditions of 
Corollary 2: 

Vs>Ai 3 t'>Ti(s,xo(Ai)) : t) > /3o(xo(Ai), t - Ti(s, xo(Ai)) V t > i' 

so the switching condition to change from mode I: x = fi(a;,d) to mode 0: a; = /o(aJ,d), is that \y{t)\ G 
[0, Ai). It is important to remark that these conditions are only sufficient and the constant Ai is, in general, 
not the threshold that defines the switching rule which leads to locally optimal performance, in the sense of 
fastest convergence. However, the fact that there exists a finite Ai satisfying the conditions of Corollary 2 
for this particular situation, means that switched control leads, if not in general to optimal performance, to 
an improvement. 
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To see farther, let us consider the following optimization problem: for the case of two systems, let 
1 3> e > be a given tolerance level with respect to which to measure convergence i.e., let us assume that 
it is required to find the shortest time for the output trajectories to satisfy \y{t,XQ,0)\ < e. Then, 

Ae(s) = argmin{ Ti{s, A) + ro(A, e) } (27) 

provides the optimal switching threshold we are looking for. If the solution of (27) is constant i.e., if A£(s) = 
A, then we set Ai = A; if A£(s) varies with the initial state value then the choice Ti(s) = T\{s , IS.s{s)) 
ensures that the output trajectory has an optimal convergence to the level e for all output trajectories of the 
switched system (4), (5). 

It may be apparent from this discussion that the difficulty of the latter optimization problem increases 
geometrically with respect to the number of systems, N, and partitions, M. Yet, satisfactory results may 
be obtained following sensible considerations on a case-by-case basis as it is further illustrated through 
particular applications in Section V. 



B. Hysteresis supervisor 

We now trade the dwell-time condition by a hysteresis assumption. Opposite to the previous case in which 
switches may occur to any mode provided that a minimal time passes, we now assume that switches occur 
as soon as the output value (in norm) leaves a determined interval modulo a hysteresis zone to prevent 
infinite switches over finite intervals. In particular, switching may occur from mode q to modes g — 1 or 
q + 1 only. The hysteresis supervisor is defined as follows: 

to = 0, i(to) = r, r G {0, . . . , M}, \h{x{tQ))\ G [A„ A^+i) ; (28a) 

4 :=arginf{|^(a;(t))| G [A^, xe, (A^+i)) } , k {i{tj) - l,i{tj) + 1} 
t>t, 

ti+i= min jtij; (28b) 

^ ke{i{t,)-i,i{t,)+i} 

i{tj+i) = k such that \h{x{tj+i))\ E [A^, xe. (Afe+i)) C [Afc, A^+i) . (28c) 
where tj, with j = 1, 2, 3, are switching times; j is the number of the last switch and the signal i{t) has 
constant value in the so-called "hysteresis" set = |J [xe (Ag+i), A^+i). 

q=0 

Theorem 2 Let Assumptions 1 and 2 hold. Then, for the system (4) with supervisor (28), measurable 
disturbances d G -Mr™ and initial conditions to = 0, xq W^, items (i)-(iii) of Theorem 1 hold with 

\y{t,Xo,d)\<p{m&x{AM,\Hxo)\},0)+-f{S[d,0,t]) Vt > to, (29) 

instead of (6a). 

Remark 5 It is important to stress that in the previous statement Assumption 1 needs not to hold globally 
i.e., each of the systems in (4) is not required to be SIIOS for all initial states and all measurable disturbances. 
It is sufficient that each system 6q is SIIOS for initial states in the set corresponding to output values where 
the system is to be active. To illustrate this, consider Figure 1 and system 6i; for the hysteresis supervisor it 
is enough that the system 6i be SIIOS for all xq such that |/i(a;o)| G [Ao, Ai] U [A2, A3]. This is of obvious 
interest if system 61 corresponds to a plant in closed loop with a controller which guarantees stability in 
the large (i.e., in a "large" specified domain of attraction, subset of M") or even locally. In Section V-B we 
present a case-study presenting such characteristic. 

Proof of Theorem 2. From (28) i(to) = r, to = where |/i(a;o)| G [A^, A^+i), r G [0, M]. According to 
Assumption 1 system 6.^ from family (1) is forward complete so estimate 

|y(t, xo, d)| < pe,,,J\h{xo)\,t) + ^{S [d, 0, h]) 
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holds also for the trajectories of system (4), (28) over t G [to,ti). From (28) the time instant ti is the time 
instant when the output trajectory (in norm) enters in the interval [ Aj^^^), X6»,(tj) with i{ti) G 

{ ^i{to) ~ 1' ^i{to) + 1 }• Assume that ti < -\-oo (if ti is infinite, then the system is clearly forward complete 
and moreover the bound (7) holds) then from forward completeness of i; = (a;,d) and considering 
d G MiRm we obtain that there exist Xq G M4. and Dq G M+ such that the solutions of (4), (28) are defined 
over [to,ti) and, moreover, 1 1 1 1 r ti) < ^0, l|d||[j„,ti) < Dq. Hence, 

-^0 = sup |fe,,,^)(a;,d)| < 00 

X£N, \X\ <Xo, |d| <Do 

which implies that for \y(t)\ generated by system Oi(^t^^) to reach adjoining intervals { [Ag.^^^^^^i, Ae^^^^j) and 
[ Ag.jj^^j+i, Ae.jj^^j+2) } it is necessary a time proportional to the maximum "speed" Fq and ti — to units of 
time. Moreover, since Fq < 00 necessarily ti > to therefore, the solutions of system (4), (28) are defined 
and are continuous over [toi^i]- Reconsidering the initial time to be ti and (4), (28) we obtain that the 
solutions of the latter are defined for all t G [ti,t2) and satisfy 

\y{t,x{ti),ti)\ < /3e,^^J\y{h)\,t-h) + -f{S[d,h,t2]) 

and there exists finite numbers Xi G M4., Di G M+ and Fi G M+ such that 

||a^||[ti,t2) < -^1, < ^1, Pi= sup |fe.„^,(a;,d)| 

xeAf, \X\ <Xi, |d| <Di 

hence t2 > h. Repeating these arguments for arbitrary j > it is possible to prove the existence of Xj G M+, 
Dj G M+ and Fj G M+ with properties 

\\x\\[t^^t,+,) <Xj, \\d\\[t^^t^^^) < Dj, Fj= sup \fg {x,d)\, 

X&Af, \X\ <Xj,\d\ <D, 

SO tj^i > tj. This implies right continuity of the switching signal i{t) and forward completeness of the 
system (4), (28). 

The rest of the proof follows along similar lines as the proof of Theorem 1: the estimates (11) and (12) 
continue to hold verbatim and if at tf^ we have |t/(tfc)| G [Am, Aj\/+i) then (13) holds true. In view of (28) 
we have either \y{tk)\ = \y{to)\ and k = or \y{tk)\ = Am (this is because according to the definition 
of the hysteresis supervisor, there is no dwell time); the bound (29) follows. The proof of statement (ii) 
follows as in Theorem 1 by replacing Tmin with 0; finally, the proof of (iii) follows the same arguments as 
Theorem 1. ■ 

Corollary 3 Let all conditions of Theorem 2 hold and let Tg, for < g < M, be defined by (23), (24). 
Then, the time T0.5 Ai that is required for \y{t)\ to converge to the set {\y\ < 0.5 Ai} for the case d{t) = 0, 
satisfies 

r 

7b.5Ai < Tr{Ar+l) + Tr+l{Ar+l) + '^Tk{Ak^i) . 

k=0 

Proof. If /i(|ic(0)|) G [x0^(Ar+i), Ar+i) then \y{t)\ may, in general, reach the interval [ A^+i, Ar+2) even 
in the case d(t) = 0. From (28) and Assumptions 1, 2 it follows that \y(t)\ may not increase and the 
trajectories of the system (4), (28) behave similarly as the trajectories of the system (4), (5), (23), (24). This 
confirms the corollary. ■ 

The advantage of the hysteresis supervisor (28) over the dwell-time supervisor (5) lays in the fact that 
for the latter, during dwell-time, the system's output can in general reach any interval [xe (Aq-|_i), A^+i), 
q G [0, Af] before a switch occurs. In contrast to this, in the case of the hysteresis supervisor (28) only 
adjoint intervals may be taken into account which results in a simpler analysis of the system's behavior. 
Additionally, the supervisor (28) ensures that only the system 9q may become active for output values in the 
interval [ A^, A^+i), while for the dwell-time supervisor (5) it is possible that the interval is reached during 
the dwell-time period by any (sub)system from the family (1). The latter allows to relax Assumption 1 for 
the hysteresis supervisor, as is explained in Remark 5. 
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On the other hand, as it is established in Corollaries 1 , 2 and discussed below the latter, a proper choice of 
output depending dwell-time functions Tq{-) can ensure that the system posses additional (desired) stability 
properties; for instance, the upper estimate (25) given in Corollary 1 for supervisor (5) is better than the 
estimate provided by the conditions of Corollary 3 for supervisor (28)). Additionally for appropriate (large) 
values of dwell-time functions Tq the switched system may admit an SIIOS like estimate -cf. [42]. Another 
shortage of supervisor (28) with respect to (5) is that the hysteresis-based supervisor allows, in general, for 
regimes with "high" rate of switches. 

Finally, it is worth pointing out that an important difference between the hysteresis supervisor proposed 
above and others published in the literature -cf. [25], is that the hysteresis and partition properties are used 
without restrictions on dwell time values. 

V. Applications 

A. Hybrid observer design for Lur'e systems 

The problem of observer design for nonlinear dynamical systems has been one of the centers of attention 
of the control community during decades and still has not been completely solved for nonlinear systems 
-cf. [44], [45]. Here, we address the observer design problem for Lur'e-type systems via the supervisory 
control approach proposed above. Consider the system 

i; = A(y)a; + (/.(y) + Bd, y = Cx (30) 

where x G i?", y G BP, d G i2™ are state, measured output and disturbance vectors correspondingly. This 
problem is well investigated e.g., it can be solved via linear systems observers design theory. Indeed, after 
[46], assume that there exists a continuous matrix function K : i?^ — > R^^p such that for some positive 
definite matrix P with dimension n x n for any y ^ BP have 

G(y)^P + PG(y) < -aP, G(y) = A(y) - K(y)C (31) 

for some a > 0. A globally exponentially convergent observer for system (30) takes form 

z = A{y)z + cl){y)+J^{y)[y-Czl (32) 

where z G BP is a vector of estimates of x. The observation error e = x — z possesses the following 
dynamics: 

e = G(y)e + Bd, 

whose origin is exponentially stable for d = 0, due to (31). In the terminology of the paper, the system 
(30), (32) is SIIOS with respect to output e and input d ^ 0. 

An unclear yet natural question is that of assigning the gain function K to ensure the best (trading off, 
speed and peaking) convergence of error e to zero. Assume that the matrix inequality (31) may be solved 
for any < amin < o < Omax then it is possible to design a gain function K to provide "slow" (amin) or 
"fast" (omax) speed of error convergence to zero. Choosing K for the "fastest" solution may, in general, 
result in large overshoots (peaking). Such negative behavior is significant for instance in the context of 
observer-based designed synchronization -cf. [47] or when the state variables z in (32), as well as inputs 
^{y) [y ~ C z], are physically meaningful. On the other hand, for small deviations i.e. , when the error e 
stays under some admissible bounds, the application of "fastest" gain K is also undesirable. Thus, for small 
and large amplitudes of observation error e it is desirable to apply a gain function K close to aniin> while 
on average amplitudes of e the "fastest" gain K may provide better performance. This leads to a natural 
switching design between "locally optimal" controls. 

To deal with this problem the proposed approach can be applied. If the system (30) is forward complete 
for any d G Mr^, then for any triplet (Kj(y), Pj, Oj), I < i < N such that (31) is satisfied the system 
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(30), (32) is forward complete and SIIOS with respect to the output e and the input d and 



2 Amax(Pi) ^ g-0.25 Oi r 



'^min (Pi 



7(s) = 2pi?max(B) s, p= max 

l<j<A' 



^max(Pfi) 
^Amin(Pi) 



where Amin(P)> Amax(P) are minimum and maximum eigen-values of matrix P and T?max(P) is the 
maximum singular value of P. Therefore, Assumption 1 holds. Assumption 2 is also satisfied for a suitable 
partition of the interval [0, + oo): in this case the supervisor (5) provides a desired switching between 
observers with different gains, which allows to avoid the peaking phenomenon while decreasing "control 
energy" spent on observation (respectively in synchronization). Next, we consider a concrete brief example; 
for simplicity, we set Tj(s) = Tmin > that is, we use a constant dwell-time. 

Example 2 Let the system (30) be a lossless pendulum i.e., 



xi =X2; y = xi; 

±2 = —J^ sin(xi) + d , 



where lij G i? is a known parameter. 




pendulum phase curve 



20 



40 




10 



1 1 

/ hybrid 


FM median 


slow 


fast 1 


1 



Xl, Zl 



4 t [sec] 



4 t [sec] 



(a) (b) (c) 

Fig. 2. Lur'e system: (a) pendulum and observer trajectories; (b) position estimation errors; (c) velocity estimation errors 



The observer (32) takes the form 

ii = Z2 + ki{xi - Zl); 

Z2 = -oj'^sm{xi) + k2{xi - Zl) , 

where K = [ki k2]'^, ki > 0, k2 > 0. Choosing ki = Xi + A2 and k2 = Xi X2 it is possible to assign any 
real negative eigen-values — Ai, — A2 to the poles of the closed-loop system. We compute three choices of 
observer gains for the system {N = 3) as Ksiow = [2 1]^ (Ai = A2 = 1); '^median = [6 9]"^ (Ai = A2 = 3); 
K/ast = [10 25]"^^ (Ai = A2 = 5), which provides different speeds of observation error convergence (for 
u = 1) -cf. Figs. 2b and 2c. The hybrid observer is constructed as follows: let the following partition of 
M+ be given (M = 3): 

Ao = 0, Ai = 0.1, A2 = 2, A3 = 5, A4 = + 00; 

with the corresponding operation modes 6*0 = median, 61 = fast, 62 = median, 6^ = slow and constant 
dwell-time Ti{s) = Tmin := O.Olsec. In view of the responses depicted in Figs. 2b and 2c we identify slow, 
median and fast modes. In Fig. 2a one can see the phase curves of the pendulum {i.e., in the space X1-X2) 
and of the observer {i.e., in the space Z1-Z2). As one appreciates, setting high observer gains for increased 
speed entails large overshoots and vice-versa. Based on such observations, to avoid peaking while conserving 
a fast convergence rate, we decide to apply the slow observer for large observation errors and to use the 
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fast observer only for sufficiently small values of observation errors. Besides, in the simulation scenario 
we have added the external disturbance d{t) = 0.05 sin(0.3t); the system's response to the perturbation is 
appreciated in the loops (closed curves) that appear in Fig 2a i.e., for "large" values of the disturbance, the 
pendulum undergoes full revolutions while for small values of the disturbance, only small oscillations are 
observed. The performance improvement via the hybrid observer is clear from Fig. 2. 

B. Hybrid synchronization of Lorenz oscillators 

Consider the problem of master-slave synchronization -cf. [48] of two Lorenz systems: 

a {yi - xi) + di{t) 

xi{p- Zl) -yi + d2(t); (33) 
xiyi - I3zi +d3{t) 

a{y2- X2) 

X2{p- Z2) - y2 + ui (34) 

X2 y2- PZ2 +U2, 

where {xi,yi,Zi) G R^, i = 1,2 are the state variables of systems (33), (34); Uj G i?, i = 1,2 are 
control inputs to system (34), u = ( ui ^2 ); dj G i?, i = 1, 2, 3 are disturbances on the system (33) hence 
d = {di d2 ds) and d G Rrs ; the parameters o" > 0, p > and /? > are identical for both systems. 
System (33) plays a role of "master" and system (34) is the "slave", it is assumed that with given disturbing 
inputs all solutions of system (33) are bounded for all t > 0. We also assume that state variables {xi,yi, Zi), 
i = 1, 2 of the master, (33), and the slave, (34), are available from measurements. It is a well known fact 
that for a = W, P = 8/3 and p = 28 with dj = 0, i = 1,2,3 the system (33) is chaotic -cf. [49]. We 
suppose that the disturbances d are small enough and they do not destroy the natural strange attractor of the 
system for the case without the disturbances. The problem of synchronization of Lorenz systems has seen a 
lot of attention during the last decades -cf. [50], [51], [52], partly because its applicability in encoding and 
secure telecommunication -cf. [53]. Briefly, in such context the goal is to generate a chaotic carrier signal 
by the transmitter (the master system) that may be decoded provided that the receiver (the slave system) 
synchronizes its dynamics (at least with respect to an output) with the master. 

While the synchronization problem may be solved under different control theory viewpoints as observer 
design -cf. [47] or tracking control -cf. [54] we wish to illustrate here how hybrid control may be used to 
achieve controlled synchronization while minimizing the input energy. To that end, we define the synchro- 
nization errors as ei = xi — X2, 62 = yi — 2/2, 63 = zi — Z2 which obey the dynamics: 

ei = a{e2 - ei) + di{t) 

62 = eip - XlZi + X2Z2 - 62 + d2{t) - Ul (35) 

63 = xiyi - X2y2 - (3e3 + dsit) - U2. 

We design control inputs that ensure input-to-state stability for system (35); in this case, systems (33), (34) 
are SIIOS with respect to the synchronization error e = (ei, 62, 63)^ as output and the input d. In particular, 
we use the hysteresis supervisor (28). 

To comply with Assumption 1 we first design two controllers that achieve the synchronization goal: 
the first is a cancellation control law and the second is based on the linearization around the origin. The 
cancellation law is given by 

ui = (^p + a - 2^^ {I - \ )a - Z2^ ei, U2 = y2ei, < A < 1, (36) 

for which system (35) has Lyapunov function V{e) = 0.5 e^e and 

V<-Xael- (v/(l-A)fTei - 62)^-/33 ei + e^d. 
Clearly the system with such Lyapunov function globally satisfies the Assumption 1. 
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The second control is local in the sense that it may be applied only if |e| < e, where e > is given: 

Ml = a ei, U2 = 0, a > 0. (37) 

With control (37) the system (35) can be rewritten as follows 






e = Ae + D(t), A 



—a a 
p — a —1 











di{t) 

d2{t)-Xi{t) Zi{t) + X2it) Z2it) 

d3it) + xi{t) yi{t)-X2{t)y2{t) 



where | |D| | < | |d| | +4 and D G Mrs is a "new" bounded disturbance. The system is linear time invariant 
and the matrix A is Hurwitz for a > 27. The system with control (37) satisfies Assumption 1 locally: on 
< e}. If e < 1 the system is input-to-state stable with input D. 

A third control input is considered which consists in applying. . . no input i.e., 



Ui =0, U2 = 0. 



(38) 



The motivation for this control law is that both, master and slave, systems have bounded solutions. For 
the former i.e., (33) it holds by assumption while for the slave system (34) we know that all trajectories 
approach a strange attractor which is strictly contained in a compact. 

We can now rewrite the system (35) as 



e = A'e + D'(t), A' 



— a 


a 








-1 














D'(t) 



diit) 

d2{t) - Xi{t) Zi{t) + X2{t) Z2{t) 
+p[xi{t)-X2{t)] 

d'i{t)+xi{t)yi{t)-X2{t)y2{t) 



with disturbance D' G M]r3. Matrix A' is Hurwitz and the system satisfies Assumption 1 globally. 

Thus, we have three systems in (4) which are composed of systems (33), (34) with controls (36)-(38). Let 
/ = { 1, 2, 3} and consider the partition of defined in Assumption 2 with M = 3 and 6q = 3, 0i = 2, 
^2 = 1 and ^3 = 3. The practical motivation for such a partition is the following: we wish to diminish the 
amplitude of control energy amplitude therefore, permanent application of the cancellation controller (36) is 
not desirable; instead, we would like to switch to the linear local control (37) for relatively median values 
of the output errors (hence A2 = e). When synchronization errors are considerably small we may afford to 
switch off the control action i.e., "control" (38) is active for output errors smaller than Ai ^ e < 1. On the 
other hand, in the case when differences in the state trajectories are "large" the use of the cancellation control 
law (36) may lead to large overshoots in control effort which is obviously undesirable since, in particular, 
may cause actuator saturation. Hence, in this situation we also use control (38) and let the state trajectories 
of the forced and unforced Lorenz models (33) and (34) converge to their common strange attractor without 
control effort. In view of the latter, A3 is chosen proportional to the diameter of a sphere strictly containing 
the strange attractor which may be computed numerically^. 

Example 3 Let us consider a numerical example. Let o" = 10, /3 = 8/3 and /? = 28; A = 0.1 and 
a = 28; Ai = 0.1, A2 = 1 and A3 = 5. For illustration, we have performed simulations applying 
the hysteresis supervisor and the controls (36)-(38) individually. For the sake of comparison we use the 
following performance functional: 

T T T 

Je = T-^ j Ht)\^dt, Ja = 10 T-^ j \e{t)\'^dt, Ju = T-^ j \u{t)\'^dt, 

.9T 

where T > defines the length of simulations' windows. The functional Jg describes the overall quality 
of synchronization in terms of the integral square error (ISE), the functional Ja gives the ISE for the last 



We are not aware of any work computing analytically the "size" of the Lorenz attractor. 
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tenth part of the simulation window, the functional Ju estimates the input control energy. We use two sets 
of initial conditions (everywhere initial time is zero): 

xi(0)=0.1, yi(0) = zi(0) = 0; y2(0) = l, ) = ) = -1; (39) 

2;i(0) = 0.1, yi(0) =zi(0) = 0; ^2(0) = 10, ) = ) = -10, (40) 

which correspond to "small" and "large" initial deviations of the systems (33) and (34), and disturbances 

= 5 sin(0.5O, d2(t) = -5 cos( 0.1 1 ), d3(t) = 2.5 sin( t ). (41) 

The simulation results are presented in Table I for the case without disturbances and in Table II for the 
disturbances (41) with T = 30 sec. As it is possible to conclude from these tables the supervisory control 
ensures the best asymptotic performance minimizing the control energy. While the overall transient is worse 
for large initial conditions (due to the application of control (38) which by itself does not solve the problem 
of synchronization) the "asymptotic" quality of synchronization expressed by the index Ja is comparable 
with those provided by controls (36), (37) but with much smaller control effort, again, due to the use of 
controller (38). 





value of Je 


value of Ja 


value of Ju 


i.e. (39) 


i.e. (40) 


i.e. (39) 


i.e. (40) 


i.e. (39) 


i.e. (40) 


Control (36) 


0.019 


1.800 








5.009 


649.191 


Control (37) 


0.034 


2.578 








5.684 


346.089 


Control (38) 


541.324 


365.279 


518.694 


350.101 








Supervisory eontrol 


0.034 


141.329 


0.009 


0.009 


5.436 


3.548 



TABLE I 

Values of performance functionals for the case without disturbances 





value of Je 


value of Ja 


value of Ju 


i.e. (39) 


i.e. (40) 


i.e. (39) 


i.e. (40) 


i.e. (39) 


i.e. (40) 


Control (36) 


0.570 


2.526 


0.864 


0.864 


71.065 


675.778 


Control (37) 


0.597 


5.181 


0.335 


0.335 


50.030 


602.633 


Control (38) 


256.362 


428.067 


474.852 


375.061 








Supervisory eontrol 


0.485 


34.445 


0.317 


0.317 


45.355 


39.365 



TABLE II 

Values of performance functionals for the case with disturbances (41) 



For completeness, some plots are shown in Figures 3 and 4. The former presents the simulation results 
for the case with disturbances given by (41) and initial conditions given by (40) under control (36). Figure 4 
depicts simulation results corresponding to supervisory control under the same conditions. More particularly, 
we show in Fig. 3a and 4a the phase portraits of the master and slave systems under control (36) and 
hysteresis-based supervisory control respectively. In Figs. 3b and 3c we show the control efforts on different 
scales under control law (36); these may be compared with the supervisory control input depicted in Fig 4b. 
As it is appreciated from the plots, the application of no control i.e., (38) for the large initial errors results 
in a serious decreasing in the applied control energy Ju and maximal amplitude of |ii(t)|. The latter has the 
cost of "high" peaks appreciated in the steady-state errors -cf. Fig. 4b; yet the overall steady-state behavior 
is comparable to that under controller (36) -cf. Fig. 3c. 
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(a) Phase portraits 



0.2 0.4 0.6 0.£ 

(b) Control inputs (zoom of first sec) 



-4D 




(c) Control inputs (zoom) 

Figure 3: Lorenz systems under cancellation con- 
trol (36) 




10 20 J 

(b) control inputs 

Figure 4: Lorenz systems under supervisory con- 
trol 



VI. Conclusion 

We have addressed the problem of establishing input to output stabilization for a family of (switching) 
nonlinear systems. Two approaches are proposed and analyzed based on dwell-time and hysteresis super- 
visors. The utility of our approach is illustrated by designing a switching observer for Lur'e-type systems 
and the performance improvement is stressed for a particular example of a pendulum. A second illustrative 
application concerns the master-slave synchronization problem of two Lorenz systems. 
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